Alternating spin chains coupled by frustrated interactions in the multiferroic 

FeTeaOsBr 
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Density functional theory calculations performed for the multiferroic FeTe205Br reveal surpris- 
ingly strong super-superexchange pathways through Fe-O-Te-O-Fe bridges, implying that the mag- 
netic lattice must be treated as a system of alternating Fe chains coupled by frustrated interactions. 
We find firm support for this model in the magnetization and the antiferromagnetic resonance 
measurements as well as in the matching of the incommensurate magnetic vector. Its peculiar tem- 
perature dependence is explained by interchain exchangestriction being responsible for the emergent 
net electric polarization. 



Frustrated low-dimensional spin systems exhibit a 
plethora of exotic magnetic ground states, which are an 
exciting challenge both for theorists and experimental- 
ists P . Since quantum fluctuations enhanced by frustra- 
tion tend to destabilize classically ordered states, such 
systems often develop complex orders with broken inver- 
sion symmetry thus fulfilling the fundamental condition 
for multiferroicity. This can lead to strong magnetoelec- 
tric (ME) coupling , most spectacularly observed as 
reversal of the electric polarization with the magnetic 
field So far, strong ME coupling has been almost 
exclusively associated with transition metal (TM) ox- 
ides. However, unconventional magnetic ground states 
@ and multiferroicity 0, 0] have been recently found 
also in TM tellurite halides, where TM ions are coupled 
through complex exchange pathways involving tellurium 
polyhedra Q [l^l ■ Due to the structural complexity of 
this family a microscopic mechanism of the ME coupling 
has yet to be established. 

As a prominent member of the TM tellurite halide mul- 
tiferroics we refer here to FeTe2 05Br 0,0, which adopts 
a layered structure with crystal layers of [Fe40i6]^"^ 
tetramer clusters connected via Tc'^"*' ions [Fig. [Ifa)]. 
The negative Curie-Weiss temperature Tew = —98 K 
implies strong antiferromagnetic (AFM) interactions be- 
tween the Fe'^+ {S — 5/2) moments, while the system de- 
velops long-range magnetic order at a considerably lower 
temperature Tni^HK [1,0, [Hi. Only 0.5 K below, at 
7^2 = 10.5 K, the second transition to a predominantly 
amplitude modulated magnetic ground state with the in- 
commensurate (ICM) magnetic vector c[icm={^ 0.463 0) 
[§| occurs and is accompanied by a spontaneous elec- 
tric polarization pointing perpendicular to the magnetic 




FIG. 1. (Color online) (a) Crystal structure of FeTe205Br 
along the be plane where light yellow polyhedra denote 
[Fe40i6]^''^ tetramer units. For clarity only Fe and Te atoms 
are shown, (b) The DFT magnetic exchange network, where 
the strongest J2 and J4 (thick lines) form alternating spin 
chains coupled by weaker frustrated interactions (thin lines). 
Marked unit cell contains 8 magnetic Fe'^^ ions. 



moments and to q/cM- The magnetic susceptibility 
was explained by assuming dominant interactions within 
tetramers Q. However, the tetramer model cannot ex- 
plain the observed ICM magnetic structure, essential for 
the ME effect in this system. The complexity of the 
TM tellurite halides arises from the exchange pathways 
that may involve several atom s. This makes their deter- 
mination much less intuitive 



12Ml4| as opposed to the 



TM oxides, where dominant TM-O-TM exchange inter- 
actions can be intuitively guessed from the bonding an- 
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TABLE I. Exchange, Ji, and easy-plane anisotropy, D, pa- 
rameters given in units of Kelvins and of J2 as calculated by 
DFT calculations for GGA-ft/, C/=0, 3eV and 5eV. The last 
column correspond to the parameters obtained from magnetic 
ground state minimization and used for the AFMR simulation 
(AFMR + GS min). 
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gles [15|, |16[ . Therefore, additional experimental and the- 
oretical investigations are needed to determine the mag- 
netic exchange network, which should unveil the micro- 
scopic origin of the complicated magnetic order as well 
as its coupling to the electric polarization. 

In this letter, we report on density functional the- 
ory (DFT) calculations for FeTe205Br, which surpris- 
ingly reveal that long Fe-O-Te-O-Fe super-superexchange 
bridges are stronger than some shorter and more direct 
Fc-O-Fc ones. Consequently, a quasi-one-dimensional 
model of alternating AFM Fe chains weakly coupled by 
frustrated two-dimensional interactions has been derived 
and verified by experiments. In this picture, the ICM 
magnetic structure appears naturally. Considering the 
fact that the ICM component of the magnetic vector 
along the crystallographic 5-axis, shows a similar 

temperature dependence as the electric polarization [llj , 
we conclude that the microscopic origin of the ME ef- 
fect is the magnetostriction of the interchain Fe-O-Te-O- 
Fe exchange pathways. These involve Te^"*" ions hosting 
easily polarizable lone-pair electrons. 

Density Functional Theory.- The underlying Heisen- 
berg Hamiltonian parameters for FeTe2 0,5Br (Tablc|l| 
were determined by total energy calculations. A full po- 
tential local orbital basis set [17| and generalized gradi- 
ent approximation (GGA) as well as GGA+U functionals 
were used. The Hubbard parameter U was chosen as 3 
and 5 cV to take into account the intra-atomic Coulomb 
interactions. The obtained network of exchange interac- 
tions (Fig. [T]) was checked for consistency by N*''-order 
muffin-tin orbital downfolding [3] calculations. 

The six considered exchange couplings are all anti- 
ferromagnetic. As expected, the intracluster (within a 
tetramer) Fe-O-Fe exchange J2 is the strongest. Strik- 
ingly, the second strongest exchange is not the remain- 
ing intracluster Fe-O-Fe exchange Ji, but rather the in- 
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FIG. 2. (Color online) (a) Measured magnetic susceptibility 
(circles) and comparison to the quantum Monte Carlo (solid 
line) and the classical (dashed line) calculations for alter- 
nating 5 = 5/2 chain with J4/J2=0.6 and J2 = 16K. Inset: 
High-field magnetization measured at T = 1.5K (circles) and 
the simulation to the model described in the text (solid line). 
The dashed line indicates the position of Bsf- (b) Tempera- 
ture evolution of ESR spectra at 405 GHz for -B j | b. Arrows 
indicate antiferromagnetic resonance modes emerging below 
the Neel transition temperature. 



tercluster super-superexchange interaction J4, mediated 
through the Fe-O-Te-O-Fe bridges. This result discloses 
that in tellurite halides the super-superexchange interac- 
tions through long bridges involving Te^"*" ions are impor- 
tant and should not be neglected on the basis of simpli- 
fied structural arguments. The two dominant couplings, 
J2 and J4, thus effectively form alternating Fe'^+ spin 
chains. The rest of the considered exchange interactions 
are weaker and provide frustrated interchain interactions 
as shown in Fig. [T] We note that the DFT results leave 
some freedom concerning the overall energy scale, but 
the ratios of J^'s are expected to be subject only to small 
errors [l^. However, the relative sizes of J^'s are some- 
what dependent on the chosen value of U , with Ji and 
J5 showing the highest sensitivity. Thus, we suspect that 
they could be modified by very small lattice distortions 
anticipated at the ME transition. 

Magnetic properties.- Quantum Monte Carlo simula- 
tions of the bulk magnetic susceptibility for alternating 
S = 5/2 chains with L = 60 sites were carried out with 
the ALPS 1.3 [i^ directed loop application 21| in the 
stochastic series expansion framework ji^]. For J2 = 16K 
and the DFT determined ratio J4/ J2 = 0.6 we find a very 



good agreement between our simulations and the exper- 
imental data for B = 0.1 T and B\\b [Fig. [^a)]. A rather 
good agreement is also achieved when classical spins are 
considered. However, since the tetramer model can also 
fit the susceptibility data Q , we need an independent ex- 
perimental proof to discriminate between the two models 
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and to fix the precise values of J's. 

Magnetization measurements performed at T = 1.5 K 
for B parallel to the ICM direction (6 axis) show a lin- 
ear increase of the magnetization with B and the exis- 
tence of a clean step at around Bs = 1 1 T [inset in 
Fig. Ilia)]. Such a behavior is reminiscent of a spin-flop- 
like process, implying the presence of considerable mag- 
netic anisotropy. 

Electron spin resonance.- Magnetic anisotropics can 
be determined by electron spin resonance (ESR) (23l - 
[26j . which in the magnetically ordered state allows a 
detection of collective low-energy magnetic excitations 
at Q = 0, ± qicM - the so-called antiferromagnetic reso- 
nance (AFMR) modes 27|. Measurements at fixed reso- 
nant frequencies for magnetic fields B\\b were performed 
at the Institute for Materials Research, Tohoku Univer- 
sity, Sendai, Japan. In the paramagnetic state a strong 
ESR signal at g = 2.005 is observed, as expected for Fe'^+ 
{S = 5/2) ions. At T^i the paramagnetic signal disap- 
pears and is replaced by new resonances [marked by ar- 
rows in Fig. EJb)] shifted with respect to the paramag- 
netic g = 2.005 value. 

The new resonances below T/vi are attributed to 
AFMR modes. Resonant fields for individual AFMR 
lines strongly depend on the resonance frequency 
[Fig.jSja)], allowing us to derive their frequency- field dis- 
persions in the range 95 to 405 GHz and to 30 T at 2 K 
[Fig-ISIb)]. The lowest and at the same time the most 
intense excitation mode marked by arrows in Fig. [3l^a) 
shows a zero-field gap Ai/^i? ^ 240 GHz [Fig. [3l[b)] cor- 
roborating a sizable magnetic anisotropy and in agree- 
ment with high-field magnetization measurements. With 
increasing magnetic field the gap reduces to ^100 GHz 
at BsF, where the slope is reversed and the gap increases 
again with increasing field. In addition, at least five more 
high-frequency modes were detected [Figs.[2]Jb) and[3]Ja)], 
which also dramatically change their behavior at Bsf- 

Simulations.- In order to relate the high-field magne- 
tization and the AFMR results to DFT calculations the 
first step is to determine the magnetic ground state of 
the infinite spin lattice via minimizing the classical en- 
ergy (compare secion[ST|) corresponding to the proposed 
Hcisenbcrg model [Fig. [IJb)]. Since Fe'^"'" ions are in the 
high-spin S = 5/2 state, such a classical treatment of the 
model is reasonably justified. In our approach we assume 
coplanar ordering and fixed sizes of the magnetic mo- 
ments thus ignoring their amplitude modulation estab- 
lished from neutron diffraction experiments Q. Starting 
from the parameters obtained from the GGA+U calcu- 
lations with [/ = eV, 3 eV, and 5 eV (Tabled, the 
solution immediately converges to an ICM order. For 
[/ = 3 eV the experimental q^cM — 0.463 is reached by 
fine tuning of the interchain couplings, i.e., Ji, J3, and 
Jg, up to 2 % (Table|l|. We note that using the tetramer 
model parameters, i.e., |Ji| = IJ2I ^ IJ3I » IJ4I, [Js], 
1^6 1 0: '^6 were unable to generate the ICM magnetic or- 
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FIG. 3. (Color online) (a) Field dependence of AFMR spectra 
at 1.5 K for B\\b and various resonant frequencies between 
95 and 405 GHz. (b) Frequency-field dependence for intense 
(thick circles) and weak (thin circles) AFMR modes. Solid 
lines are simulations to the model described in the text. The 
dashed line indicates Bsf. 



der, which clearly demonstrates that the tetramer model 
is inappropriate. 

Due to the excellent agreement with <?fpj\/ we now 
use the optimized parameters (TablelH fourth column) 
to calculate the AFMR modes. For computing the equa- 
tions of motion [23| we define a finite magnetic lat- 
tice comprising 7 unit cells coupled along the crystal- 
lographic b axis, i.e., mimicking the experimentally ob- 
served gfp^j = 0.463 ~ 3/7. We assume 5 /is for the 
size of Fe'^"^ moments and introduce an easy plane 
anisotropy, DS^ . Here D is the magnitude of single- ion 
anisotropy. To fix the direction of the hard axis along 
z = (0.31, 0, 0.95) we take into the account the orienta- 
tion of the ordered magnetic moments Q and a two-fold 
screw axis, which coincides with the b axis. The only free 
parameters left in these calculations are thus D and the 
strength of the J2 interaction (all other interactions are 
scaled appropriately). 

The strongest AFMR mode is best described [Fig.[31Jb)] 
by J2 = 19.0K and i:) = 0.18K (TableU]). We stress that 
the theoretical curve nicely reproduces the softening of 
this mode up to Bs f and its reversed character at higher 
fields. At the same time, the dispersions of all other in- 
tense higher-frequency AFMR modes are in convincing 
agreement with the experiment [Fig. El^b)]. Moreover, 
the same set of parameters perfectly simulates the mag- 
netization response to the applied magnetic field includ- 
ing the magnetization step at Bsf [inset in Fig.[2][a)]. To 
conclude this part, we emphasize a remarkable agreement 
between experiments and ratios of J's obtained by the 
DFT calculations for GGA+U with C/ = 3eV achieved 
only by scaling slightly the exchange interactions. 

Discussion.- The main result of this study is the 
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finding that FeTe205Br has to be treated as a system 
of alternating 5* = 5/2 chains weakly coupled by frus- 
trated interactions, in contrast to the previously pro- 
posed tetramer model. The knowledge of the appropriate 
spin model allows us to investigate the ME mechanism 
from a microscopic perspective. In order to identify the 
exchange pathway responsible for the ME effect we first 
recall that q^c^ is temperature independent in the high- 
temperature ICM phase (between Tati and 2^2), while 
in the low-temperature ICM phase (below Tjv2) it scales 
similarly to the electric polarization, i.e., it behaves as 
(Tjv2 — T)'^-^^ To address this important point we 

return to the minimization of the classical energy for 
the infinite lattice and calculate how g^f^^,^ is affected 
by small changes of different J^'s. Since the ICM modu- 
lation is perpendicular to the alternating chains, it suf- 
fices to vary only the interchain exchange couplings Ji 
(i=l,3,5,6) with respect to the values that reproduce the 
AFMR modes at T = 2K. In order to reproduce the de- 
crease of q^cM observed below T/v2, Ji has to be reduced, 
while J3, J5 and Jg have to be increased (Fig. EJ. 

Given that Ji and J3 are superexchange interactions 
via Fe-O-Fe bridges, they are expected to behave ac- 

If Ji 
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cording to Goodenough-Kanamori rules [If 
is responsible for the ME effect, then the Fe-O-Fe angle 
should be reduced or the Fe-0 distance increased in or- 
der to weaken Ji. No significant structural changes for 
this bond have been observed in our high-resolution syn- 
chrotron x-ray diffraction experiments [8| (compare also 
section rS3|) thus ruling out this possibility. An analogous 
conclusion can be derived also for J3. We are thus left 
with the two remaining candidates - J5 and Jq super- 
superexchange interactions - which both involve Te^"*" 
ions. Since the only sizable change at the ME transition 
corresponds to the shortening of the Fe2-Te3 distance, in- 
volved in the J5 pathway, we assign the microscopic ori- 
gin of the ME coupling to the changes of this interaction. 
The key to understanding the ME effect is hidden in the 
low symmetry of the long Fe-O-Te-O-Fe bridges, which 
allows emergent magnetic order to provoke a softening of 
an appropriate phonon mode and hence a net electric po- 
larization. We stress that this is possible even if the spin- 
orbit coupling is very small. This resembles Ni3V208, 
where the magnetostriction mechanism is proposed to be 
responsible for the ME couping [28|. Finally, we point 
out that tellurite halide systems seem particularly in- 
clined to such effects, since their long low-symmetry ex- 
change bridges involve easily polarizable Te^'*' lone-pair 
electrons, which enable that already minimal changes in 
the strength of the exchange interactions are sufficient to 
produce a measurable macroscopic effect. 

Conclusions." We have shown that the magnetic ex- 
change network of the FeTe205Br system should be re- 
garded as a system of alternating Fe chains with weaker 
frustrated interchain interactions. Frustration is found to 
be responsible for the observed low-symmetry ICM mag- 
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FIG. 4. (Color online) Dependence of q^(^]^i on interchain 
exchange coupling constants ratios, Ji/Jf for i = l,3,6, cal- 
culated by direct minimization of the classical energy, where 
jf are exchange parameters used in AFMR simulations. In- 
set: the same dependence calculated for J5. 



netic ordering, essential for the establishment of the mul- 
tiferroic phase. Additionally, we show that in FeTe205Br 
the ME effect on the microscopic level originates from the 
magnetostriction of interchain J5 Fe-O-Te-O-Fe exchange 
pathway, where the net electric polarization comes from 
the Te^"'" lone-pair electrons. Finally, our findings clearly 
demonstrate that in tellurite halides one needs to be ex- 
tremely cautious when making assumptions about the 
exchange network, as arguments based solely on the 
crystal structure may be very misleading. Long super- 
supercxchangc bridges involving Tc^"*" ions can lead to 
surprisingly strong interactions. In addition, their struc- 
tural volatility makes the ME effect possible. 
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with the x-ray experiment and data analysis. This work 
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dation project 200021-129899, the DEC (SFB/TR 49 and 
HO 2325/4-2), and by the Helmholtz Association through 
HA216/EMMI, and the Slovenian research agency (Jl- 
2118). 



SUPPLEMENTARY INFORMATION 
SI. Magnetic ground state calculations 

We consider initially the simple Heisenberg model, 

= ^ JijSi ■ Sj, where Jij= Ji,..., Jg. (1) 

Due to the hard axis anisotropy, we will consider in 
what follows coplanar moments. In this case the classi- 
cal energy can be written as a function of the angle 9i 
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FIG. SI. (Color online) Measured magnetic susceptibility 
(circles), classical (dashed line) and quantum (solid line) 
Monte Carlo results for an alternating S = 5/2 chain with 
60 spins, J4/J2 = 0.6, and J2 = 16 K. Here we also include 
results for a two-dimensional 64 x 64 lattice of classical spins 
(dotted line), using the exchange ratios determined from the 
ground state minimization and AFMR calculations, but with 
D = and scaled to J2 = 12 K. Statistical errors of the 
Monte Carlo sampling procedure are well below the width of 
the lines. 



corresponding to the orientation of the spin i with re- 
spect to the predetermined direction within the confined 
plane 



ci 



(2) 



The classical energy thus depends only on the differences 
of orientations from spin to spin. 

For a spiral order, Eci can be written as a function 
of 8 angles 6i and the twist angles and defined as 
follows: 



a 



^i,nry~\-mrz 



(3) 



where n and m are the cell indices. 

We used the conjugate gradient method to minimize 
the energy and to find the magnetic structure of the clas- 
sical ground state. 



S2. Magnetic susceptibility calculations 

The magnetic susceptibility was computed by classical 
and quantum Monte Carlo simulations of the Heisenberg 
model dl]) and compared to our single crystal x(T) mea- 
surements for B = 0.1T and B\\b. Results for the alter- 
nating spin chain model were already presented in the 
main text; here we also present results for magnetic sus- 
ceptibility of a two-dimensional lattice of classical spins 



with the exchange ratios determined from the ground 
state minimization and AFMR calculations (Fig. ISip. In 
order to explain the observed high-temperature behav- 
ior within the two-dimensional model, we need to scale 
to J2 = 12 K, while for simplicity we set D = 0. All 
numerical results are for a fixed system size and periodic 
boundary conditions, but we have checked that finite-size 
effects are negligible except possibly at very low tempera- 
tures. Theoretical results were converted to experimental 
units assuming a spectroscopic g- factor g = 2. 

Inspired by the scaling of the high-temperature behav- 
ior with spin quantum number S, we used a phenomeno- 
logical scaling 



XS,phcn(r) ^ XcliSiS + 1) T) 



(4) 



to map the classical result Xci to the quantum result for 
spin S. As is illustrated by the results for the alternating 
chain, this works well at high temperatures, but misses 
a suppression of x by quantum fiuctuations at low tem- 
peratures. 

In the two-dimensional model, interchain coupling is 
frustrated. Since this gives rise to a sign problem in 
the quantum Monte Carlo approach, we have to rely 
exclusively on the classical Monte Carlo simulations for 
this case. The two-dimensional model fits the experi- 
mental magnetic susceptibility at high temperatures at 
least as well as the chain model, but the agreement is 
worse around and below the maximum of x{T) (Fig. lSlj ). 
However, we have to keep in mind that in analogy to 
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FIG. S2. Refinement results for the 0-Fe distances obtained 
from synchrotron x-ray diffraction, corresponding to (a) Ji 
and (b) J3 exchange pathways. 
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the one-dimensional model, quantum fluctuations are ex- 
pected to reduce x(^) lo"^ temperatures also in two 
dimensions. Furthermore, the frustrated nature of the 
interchain coupling might enhance this reduction. In 
any case, we attribute the fact that a one-dimensional 
model yields a good effective description to the frus- 
tration in the interchain coupling. Finally, we note 
that the two-dimensional model exhibits signatures of 
an ordering transition at T « 9.5 K. This is remark- 
ably close to the experimentally observed ordering tran- 
sitions, although the Mermin- Wagner theorem forbids a 
true finite-temperature ordering transition in a strictly 
two-dimensional Heisenberg model. 

In conclusion, the fact that one-dimensional, two- 
dimensional and tetramer [§] models provide good fits to 
the magnetic susceptibility of FeTe205Br demonstrates 
that further data is needed to clarify the nature of the 
microscopic exchange network. 



S3. Synchrotron x-ray diffraction 

Low temperature single crystal X-ray synchrotron 
diffraction measurements were performed at the BMOlA 
Swiss-Norwegian Beamline of ESRF (Grenoble, France). 
Data sets (typically 780 reflections per temperature 
point) were collected in the temperature range 4.5 K 
to 35 K at a wavelength of 0.64 A, using a closed-cycle 
He cryostat mounted on a six-circle kappa diffractometer 
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KUMA. Here we show refinement results for interatomic 
distances corresponding to interchain exchange pathways 
(Figs. [S2l [S3| . The only significant changes occurring 
below T/V2 = 10.5K, involve Tea, i-G-i the J5 exchange 
pathway. 



FIG. S3. Refinement results for the Fe-Te distances obtained 
from synchrotron x-ray diffraction, corresponding to (a) Jg 
and (b) J5 exchange pathways. 
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